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Abstract 

We present here three new solutions of Brans-Dicke theory for a stationary 
geometry with cylindrical symmetry in the presence of matter in rigid rotation 
with T^n y£ 0. All the solutions have eternal closed timelike curves in some 
region of spacetime the size of which depends on uj. Moreover, two of them 
do not go over a solution of general relativity in the limit uj — > oo. 
PACS number(s): 04.50+h, 04.20Jb 
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General Relativity (GR) has survived many experimental tests during its rather long 
existence. However, there remain some conceptual difficulties (mostly related to Quantum 
Gravity) that have led to the construction of several alternative theories (for a review, see 
]1|]). Brans-Dicke theory (BDT) |2|] (along with its generalization, namely the family of 
scalar-tensor theories ||) seems to be the most promising of these because of its excellent 
agreement with experiment and its intimate relationship with the low energy limit of string 
theory [[|. To be more precise, BDT is consistent with solar-system experiments if the 
coupling parameter satisfies the inequality |u;| > 500 f|. As a consequence, a great deal 
of effort has been devoted to explore the whole space of solutions of BDT. We will be 
concerned here with solutions that display causal anomalies in the form of closed timelike- 
curves (CTCs). It has been known for a long time that CTCs are a distinctive feature of 
geometries in which the light cones may be tilt in an appropriate fashion. The most widely 
known geometries of these type are the so-called wormholes || which are exact solutions of 
a given theory of gravitation with non-trivial topology (see for wormholes in GR, and || 
for the case of BDT). They do not exhibit CTCs ab initio; some (rather artificial) procedures 
to generate them are explained in ||. On the other hand, the desired tilting is naturally 
achieved in rotating space-times such as the spinning black hole of Kerr [O, the infinite dust 



cylinder of van Stockum |TT|, the universe of Godel [0], and their generalizations Hl3| . In 



this paper we present solutions to BDT for a cylindrically symmetric spacetime with rigidly 
rotating matter modeled by a stress-energy tensor with nonzero tangential stresses. We shall 
see that the solutions have eternal closed timelike curves in a certain region of spacetime, 
the boundary of which depends on the value of u. We also analyze the singularity structure 
of the solutions. Finally, we study the u — > oo limit, and we find that two of our solutions 
do not go over solutions of GR. 

If for simplicity units are chosen which make the velocity of light equal to unity, the field 
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equations of BDT are [] 



where G^ v is the Einstein tensor, T Mi , is the energy-momentum tensor of matter and all 
nongravitational fields, and T = T 7 ^. In a stationary spacetime the metric and the scalar 
field can be chosen in such a way that g^o = t o = 0, and g^ = 0. Further simplification is 
possible if we assume that the geometry has cylindrical symmetry. The line element can be 
written in this follows: 

ds 2 = dt 2 - e 2X(T \dr 2 + dz 2 ) - l(r)d0 2 + 2m{r)dBdt. (3) 

We shall take for the stress-energy tensor the expression, 

Tp, = pu^ + (4) 

with p the density of energy, and = 5^ the four velocity in the comoving system. is the 
anisotropic pressure tensor, which is symmetric, trace-free, and orthogonal to the comoving 
observer. We will restrict ourselves here to the case in which = diag(0, —a, a, 0) 0. 
With the foregoing assumptions, the nonvanishing field equations are 



x In what follows, Greek indices run from to 3 while Latin indices from 1 to 3. Semicolons 
denote covariant derivative with respect to the metric g^, primes derivatives with respect to r, 
and the commas mean partial derivatives with respect to the coordinate x^. As usual, g = det g^ v 
and 5"^ is the Kronecker delta. Sub-indices with hat refer to the proper reference frame, in which 
g^t, = g^, being g^ v the metric of Minkowski spacetime. 

2 Let us remark that a solution for the same geometry in the case of dust in BDT was found by 



Bandyopadyhay [14], but it exhibits the undesirable feature of a singularity at a finite radial proper 



distance from the origin. 
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D 2 = / + m 2 . 



(5h) 



To solve this system, first note that Eq. (pel) readily gives the first integral 



ml b 
2D ~ 6 



(6) 



where b is an integration constant. From Eqs. (|5a| ) and ( |5gD we obtain at once that 

2D6 2 



(2 + ^)v^D^ 



(7) 



whereas from the definition of □</> one has 

2 6 2 + (2 + cj)00" 



.D = exp | — J 



dr 



(8) 



(2 + cj)00' 

Eq. QSclD in conjunction with Eqs. (phi), (Q) and (§) determine a differential equation for 
the BD scalar. We shall propose three different types of solutions. In the case of Type I 
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solutions we assume that the scalar field is given by a polynomical expression in r; in Type 
II, by a constant times a trigonometric function, and in Type III we assume that <fi is given 
by a constant times an hyperbolic function. After an appropriate fixing of the constants to 
recover GR whenever it is possible, we obtain the expressions of 4>(r), l(r) and m(r) listed in 
Table I, along with the range of values of u in which each solution is valid. Finally, adding 
Eqs. (fjeT) and (||) we obtain an equation for A, 

= y/=^ + !) □</> _ _ 1 r _ D»_ 

D <p 2 (j) 2 2 cf> AD 2 2D U 

Integrating this we obtain the expressions listed in Table II. 

In order to analyze the singularity structure of these spacetimes, we next study the 
behaviour of certain scalars built from the Riemann tensor. It is easier to compute this 
tensor in the orthonormal frame. Let us define the differential forms 

6° = -m d6 + dt 6 1 = e x dr 6 2 = e x dz 9 3 = D dQ (10) 

with the corresponding basis vectors 

4 = b\ e\ = e~ x 5\ e\ = e~ x 5\ e\ = D~ l 5\ + mD' 1 8% (11) 

The relationship between the Riemann tensor in the two frames is the usual one: 

In the following we shall see that the first two solutions of Table I are unphysical. From 
the expression of for the polynomial case, it can be seen that there is a singularity 

at the origin unless T = 0. This implies that a = 0, and so the solution reduces to 
the dust case previously analyzed in It is worth recalling that in this case there 

is a singularity which occurs at a finite proper distance from the axis of symmetry, since 



'Actually, one still has to require that spacetime be Euclidean near the origin. This fixes C to 



zero. 
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J °° exp(X)dr is convergent. Solutions of Type II has instead naked singularities for r = mr/2, 
with n G N '. Consequently, from now on, we shall concentrate only on Type III solutions 
(valid for u G (—00, —2)) for which the components of the Riemann tensor are well behaved 
throughout the entire spacetime. 

The result of the calculation of the kinematical quantities associated with the fluid shows 
that the expansion scalar, the acceleration vector, and the shear tensor are all zero, but there 
is a nonzero vorticity vector given (in the case of Type III solutions) by 

u)» = (o,0, -^/M/2- l(coshr)- 3 - 2 ^ e^' 2 >\o\ (12) 

Let us turn now to the analysis of the relation between cause and effect in our spacetime. 
First, it is necessary to remark that, although in the metric of Type I the coordinate 9 is 
an angular coordinate, this is not the case in Type III. However, after a trivial coordinate 
transformation the line element may be written in the form 

ds 2 = dt 2 - e 2A(r) (dr 2 + dz 2 ) - L(r) d$ 2 + 2M(r) dt d$ (13) 

where 

A(r) = (1 + u) ln(coshr) - (2 + ui) r 2 /A, 
L(r) = F 2 (r)[l + 2(2 + u) sinh 2 r] sech 2 r, 

and 

M(r) = ^2\u\ - 4 F(r) tanhr, 

where F(r) is any function that goes to as r 2 . The coordinate $ is now an angular 
coordinate. Consequently, any curve with constant t, r, and z is closed. In particular, such 
closed curves are timelike (but they are not geodesies) if r > r crit = arcsinh[l/(4 + 2a;)]. 

In order to characterize in more detail the matter that generates the geometry, we deter- 
mine next whether it violates the weak energy condition [HJ To do so, we calculate a(r) 
fromEq. (H) 

4 We stress that we are concerned here only with the matter sector of the stress-energy tensor. Note 
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i- n = v 2 (M 2 ) (coshr) 2 '"'- 1 exp J-Jf^—^ r 2 ) (14) 



Finally from (pg[), p comes out as 



4vrp = 6 (2H - 3)(coshr) 2|w| ~ 3 exp (-^ r 2 ) (15) 



'2(M 

It is straightforward to compute that p > 0,Vr if b > 0. Thus matter will satisfy WEC if 
and only if p ± a > or equivalently if r > r WEC = arcosh 



2(2|oi|-3) 
|cj|-2 

From Fig. 1, we can see that the radius of the causal region decreases fast with uj, and 
that there will be CTCs both in the presence of "exotic" and "non-exotic" matter for all 
values of a;. It can also be seen that the radius of the region in which matter violates WEC 
tends rapidly to a constant. 

We close with some remarks regarding the limit uj — * oo. The study of the limits of 



geometries depending on some parameter has been initiated by Geroch [16], who pointed 
out that such limit may depend on the coordinate system chosen to perform the calculations. 
More recently, Paiva et al ]Tj]] developed a coordinate-free approach (based on the character- 
ization of a given space-time by the Cartan scalars [0 ) to study these limits in GR. Romero 



and Barros found later [19| some examples of BD solutions which do not reduce to GR when 
uj — > oo. This motivated the application of the above mentioned coordinate- free method to 
BDT by Paiva and Romero f20j . The method involves a lengthy calculation (crucial to make 
any statement about the limit uj — > oo) which will be presented elsewhere. Note however 
that the BD scalar in solutions of Type II and III goes to zero in the limit uj — > oo, and 
so the effective gravitational constant G c q = G /4> diverges. Clearly our solutions do not 
go over GR, and are then counterexamples to the recent claim by Banerjee and Sen that 
the condition T ^ is a necessary and sufficient condition for BD solutions to yield the 



however that one should take into account both the BD and the matter sectors of the stress-energy 
tensor in order to relate violations of WEC to the behaviour of geodesies. 
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corresponding solutions with the same energy-momentum tensor 
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TABLES 

TABLE I. The table shows different choices of the BD scalar 0, the corresponding metric 
coefficients, and the maximum range of uj for which each solution exists. 



(f>(r) l(r) m(r) Domain of to 

l-r 2 b 2 /(4 + 2uj) (r 2 -6 2 r 4 )[l-r 2 6 2 /(4 + 2w)]- 2 br 2 [l - r 2 b 2 /(4 + 2uj)}- 1 (-00, -2) 

(&V / 2/v / 2+~w) cosr [1 - 2(2 + a;) sin 2 r] sec 2 r ^A + 2uj tanr (-2, 00) 

(b^f^2/^/T+^) coshr [1 + 2(2 + uj) sinh 2 r] sech 2 r ,/2|£j| -4 tanhr (-00, -2). 



TABLE II. The table shows the three different expressions for A(r) in correspondence with 
those given in Table I. Their behavior with the radial coordinate is also indicated. 



A(r) Typical feature of the solution 



(1 + oj) ln[2(2 + to) — b 2 r 2 ] + Tr + C Singularity at a finite proper distance. 

(1 + uj) ln(cos r) + (2 + u) r 2 /4 Naked singularities at r = mr/2. 

(1 + (J) ln(coshr) — (2 + uj) r 2 /A Well behaved throughout the spacetime. 
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FIG. 1. Variation of r cr i t (black line) and rwEC (grey line) with \lo 
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